The Cosmological Axino Problem 
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We revisit tiie cosmology of the supersymmetric QCD axion, highlighting the existence of a 
serious cosmological axino problem that is fully analogous to the gravitino problem of overclosure 
via thermal production. A general analysis implies that the QCD axino has a mass greater than 
or equal to that of the gravitino in the absence of unnatural fine-tuning or sequestering. As a 
consequence, bounds from thermal gravitino and QCD axino production are complementary in 
parameter space, and together provide a quite stringent limit on the reheating temperature after 
inflation given by Tr < 10^ - 10® GeV for an axion decay constant of fa = 10® - 10^^ GeV. Motivated 
by this result, we explore the cosmology of gravitino LSP and axino NLSP at low Tr and present 
three realistic scenarios for dark matter. 



I. INTRODUCTION 

The QCD axion is an extraordinarily elegant solution 
to the strong CP problem which is uniquely amenable to 
a variety of astrophysical and laboratory probes. Like- 
wise, weak-scale supersymmetry offers a theoretically 
motivated resolution to the gauge hierarchy problem and 
a wealth of implications for the LHC. In this paper, we 
embrace both of these theoretical proposals and explore 
the cosmology of the supersymmetric QCD axion. 

By construction, the QCD axion couples to the gluon 
with a strength inversely proportional to the axion de- 
cay constant, fa- Supersymmetry then requires a corre- 
sponding coupling of the QCD axino to the gluino. Ex- 
pressed in superspace, this interaction takes the form 



C = / d'^eAWW + h.c. 



(1) 



where the superfield containing the saxion, axion, and 
axino is defined as 



"Planck slop" operators induced by unspecified ultra- 
violet physics, they also arise irreducibly from the cal- 
culable dynamics of supergravity Q. Since these oper- 
ators transfer supersymmetry breaking effects of order 
™3/2 ~ F/mpi into the PQ sector, the axion and the 
axino acquire a mass difference of this order. In turn, 
because the axion is ensured to be massless by the non- 
linear ly realized PQ symmetry, the axino has a mass of 
order 7713/2 due to Planck-scale corrections. 

Concretely, this simple physical argument can be un- 
derstood by the existence of a higher-dimensional opera- 
tor in the effective field theory below fa^, 



,4, (A + A^riX + X^) 



rnpi 



(4) 



which exactly preserves the full PQ symmetry,^ and is 
expected from Planck-scale dynamics. Hence, the axino 
mass has a lower bound of order 



> 



(5) 



X 



A = {s + ia) I V2 + ^/2eh + Fa- 



(2) 



In this basis A shifts by an arbitrary imaginary constant 
under the nonlinearly realized Peccei-Quinn (PQ) sym- 
metry. In this basis the axino is defined to be a PQ 
singlet. While Eq. ([T]) may be the only axino coupling 
present, as in the case of the KSVZ l] axion, there can 
be additional couplings between the axino and other su- 
perpartners, as in the case of the DFSZ Q axion. 

Since the cosmology of these theories depends crucially 
on the size of the axino mass, it is of utmost importance 
to ascertain its typical value. On general grounds, one 
expects supersymmetry breaking to seep into the PQ sec- 
tor by way of higher-dimensional operators that couple 
the axion directly to the supersymmetry breaking field, 



X = x + V20r] + 9^F, 



(3) 



where F is the supersymmetry breaking scale. While 
these contributions may be present in the form of 



which provides a theoretical motivation for considering 
theories of low-scale supersymmetry breaking with grav- 
itino LSP and axino NLSP. Eq. ([5]) embodies a central 
claim of this paper: barring fine-tuning, the axino mass 
acquires an irreducible contribution of order the gravitino 
mass which is allowed by all symmetries and is generated 
by uncontrolled Planck-scale physics. 

A notable exception to this simple argument arises in 
certain scenarios in which extra dimensions are used to 
effectively sequester supersymmetric breaking from the 
PQ sector [134 ISf . Nevertheless, even in these theories, 
supersymmetry breaking is still mediated through super- 
gravity effects which typically yield the same result as 



^ Recently a similar effective theory has been considered in 
^ Explicit PQ breaking operators can very easily destabilize the 
axion solution to the strong CP problem, even if they are gen- 
erated by Planck scale dynamics See [3 4ll|l for possible 
resolutions to this difficulty. 
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in Eq. ^ [l^. In this paper we will assume that these 
more complicated dynamics are not at play. 

As we will see, the bound in Eq. ([5]) has an enormous 
effect on early universe cosmology. In particular, it im- 
plies the existence of a cosmological axino problem which 
is similar to the well-known gravitino problem [l3| but 
which occurs in a complementary region of ^713/2-^ To- 
gether, the combined bounds from axino and gravitino 
production completely exclude the possibility of a high 
Tfl, as shown in Fig.[TJ In turn, this rather unequivocally 
nullifies the viability of high-scale leptogenesis UM while 
still permitting low-scale models of soft leptogenesis [l^ 
and testable theories of asymmetric freeze-in |20l . [2lj | . 

Let us now briefly outline the remainder of this pa- 
per. In Section 2, we present the theoretical rationale for 
Eq. ^ , both in general and for a canonical supersymmet- 
ric axion model. We then go on to explore the thermal 
production of axinos in the early universe in Section 3. 
Here we take note of a novel regime in which the domi- 
nant mode of axino production arises from the decays of 
superpartners still in thermal equilibrium, i.e. freeze-in 
[22| . Freeze- in of gravitinos and more gene rally of hid- 
den sector dark matter was considered in [23j and 0,[2^. 
By combining the bounds on overclosure from gravitino 
and axino production, we can then precisely quantify the 
seriousness of the cosmological axino problem. We go on 
in Section 4 to consider the mixed cosmology of graviti- 
nos, axinos, and saxions. Applying bounds from overclo- 
sure, structure formation, and big bang nucleosynthesis 
(BBN), we determine the values of m^^2j "^a and fa for 
which the cosmological history is viable. Finally, we con- 
clude in Section 5. 



II. SUPERSYMMETRIC AXION THEORIES 

To begin, let us provide a more rigorous justification 
for the naturalness argument that implies Eq. ([S])- Ev- 
ery supersymmetric axion theory can be described by an 
ensemble of interacting fields 

= (l)^ + V20^, + e^F„ (6) 

whose dynamics induce vacuum expectation values, 
= fi- In the supersymmetric limit, the saxion, ax- 
ion and axino are linear combinations of Re (/i^, Im 4>i 
and ipi, respectively, and are massless as a consequence 
of the nonlinear ly realized PQ symmetry. The entirety 
of our discussion will occur in the language of linear 
fields, $i, so let us briefly note that the field A discussed 
in the introduction is a linear combination of the fields 
A, = Si + ^a^ + -I- &^Fa, where = fie"^^^^^.^ 



Next, let us consider the effect of higher-dimensional 
operators of the form 

which couple supersymmetry breaking to the PQ sector. 
Here is an order unity dimensionless coupling and X 
is as defined in Eq. ([3]). The parameter A is the mass 
scale of the heavy particles which couple the PQ and 
supersymmetry breaking sectors and, as discussed in the 
introduction, it is at most the Planck scale, A < mpi. 

It is clear that the only fermion mass terms gener- 
ated by this operator are a mixing between ipi and 77, 
the fermionic component of X. Morally, this can be un- 
derstood as a rather surprising mixing term between the 
axino and the goldstino. Such a mixing is strange, but 
by including a dynamical supersymmetry breaking sector 
for X and diagonalizing the fermion mass matrix in the 
presence of Eq. ([7]), one goes to mass eigenstate basis and 
expectedly finds a massless goldstino and a heavy axino 
with a mass of order 

F 

rua ^ ^, (8) 

as shown explicitly for the very simple theory described 
in Appendix A. 

One can reach the very same conclusion through a less 
technical, more physical argument. In particular, the 
right-hand side of Eq. ^ manifestly induces a non-zero 
value for the auxiliary field Fi, 

F. = -^ + ..., (9) 

which mediates supersymmetry breaking effects propor- 
tional to Fi directly into the scalar potential of the PQ 
sector. Inserting appropriate powers of the characteristic 
scale of the PQ sector, fi, we find that the mass scales 
and vacuum expectation values in the PQ sector shift by 
an amount which is of order F/A. Hence, one should 
expect a mass splitting between the axion and the axino 
of order F/A which implies Eq. ([SJ. 

This argument can be restated in another way. In par- 
ticular, apply a field redefinition 

which removes the connector interaction in Eq. ([T]) at 
leading order in 1/A at the expense of introducing X 
directly into the PQ sector dynamics. The effect of X 
will be to explicitly generate an axino mass. 



^ As in the case of the gravitino problem, this axino problem can 

be evaded if R-parity is broken. 
* Naively, the 0-^ component of 3>i contains a term quadratic in the 

axino which can produce what appears to be an axino mass in 



interactions involving only a single 4^. However, this quadratic 
term appears in the combination (-Fa; —o,i /2), and so can always 
be removed by a shift of the auxiliary field. 
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For example, consider the canonical supersymmetric 
axion theory defined by canonical Kahler potential and 
a superpotential 



W 



(11) 



where k is a dimensionless coupling, and a straightfor- 
ward calculation shows that the vacuum is stabilized at 
/1/2 = and /a = in the supersymmetric limit. ^ 
Consider the effect of the higher-dimensional operators 
in Eq. ([7]) for the symmetrical case A 1,2 = A and A3 = 0. 
After the field redefinition in Eq. PH)) . a simple calcu- 
lation shows that the vacuum is slightly shifted so that 
/s = ^XF/kA. This induces a mass for the axino equal 
to XF/2A from the superpotential, which again accords 
with Eq. (P. 

Before we continue on to cosmology, let us comment 
briefly on some of the existing literature on the mass of 
the axino. While statements are frequently made to the 
effect that the axino mass is highly model dependent, 
we disagree. The authors of 16] and [1^ computed the 
axino mass in a variety of settings and verified Eq. ([5]) 
and Eq. ([5]) in all cases. The only exceptions resulted 
from the imposition of special relations among disparate 
and unrelated parameters to yield a lighter axino. We 
conclude that in generic theories the axino mass is given 
by Eq. (|5]) and the only model dependence is the value 
of A, which has an upper bound of mpi. 



III. AXINO COSMOLOGY 

Given a proper theoretical justification of Eq. ([S]), let 
us now consider the early universe cosmology of the ax- 
ino. Like the gravitino, the axino is produced through 
thermal scattering and decay processes. As observed by 
Strumia [29i] , one can trivially compute axino production 
by translating every equation relevant to gravitino pro- 
duction by the simple replacement 



rrig V2a3 



(12) 



while shutting off all scattering and decay processes in- 
volving gauginos and scalars not accounted for in Eq. ([l} . 

In Fig. [2] and Fig. [3l we have numerically plotted 
contours of ila/i^ = 0.11 for rrig = 638 GeV and dif- 
ferent values for the axion decay constant, assuming 
the axino to be cosmologically stable. As one can see, 
for axion decay constants within the "axion window" , 
10^ GeV < fa < 10^2 GeV, the bounds on Tr are quite 
stringent. For a GUT scale axion decay constant corre- 
sponding to the "anthropic window" , the bound is weak. 



^ Our results hold irrespective of the precise dynamics which break 
supersymmetry — indeed, one obtains the correct answer even 



10' - 



when treating 
breaking. 



X as a non-dynamical spurion of supersymmetry 



10- r— 



10-' 10-5 [Q-4 [Q-3 [0-2 [Q-i [qO iQl iq2 

»i3,2[GeV] 



FIG. 1: The red contour is (13/2/1 ~ 0.11 from gravitino 
production alone, and the blue contour is flah^ ~ 0.11 for 
a^xino production alone with ma — 7713/2 and fa = 10^^ GeV. 
For larger values of ma and smaller values of fa, the blue 
contour is translated to the left. The green contour can be 
interpreted in two ways: (i) as the + ^a)h^ = 0.11 

contour for combined gravitino -I- axino co-dark matter when 
both are cosmologically stable for ma — 171^/2; (ii) as the 
^3/2/1^ ~ 0.11 contour for gravitino dark matter when the 
axino is unstable and decays sufficiently quickly to gravitinos, 
for any value of ma > m^p. Observe that a high reheating 
temperature, Tn > 3 x 10^ GeV, is unambiguously excluded 
in both cases. The superpartner spectrum is taken to be 
{m^,mi3,mg} = {100 GeV, 210 GeV, 638 GeV} and univer- 
sal scalar masses at the GUT scale equal to 500 GeV. 



Moreover, as in the case of gravitino cosmology, one dis- 
covers regions where axino production arises dominantly 
from freeze-in and scattering, corresponding to the ver- 
tical and sloped portions of the contours, respectively. 
Note that the axino contours in Fig. [5] are "fiipped" from 
the usual contours corresponding to gravitino overclo- 
sure, simply because the scattering and decay rates for 
axinos are independent of the axino mass and fixed by 

fa- 



A. Cosmological Axino Problem 

The cosmology of gravitinos and axinos has been stud- 
ied in great detail by numerous authors over many years 
[3]| . Nevertheless, to our knowledge the existence of a se- 
rious cosmological axino problem has never been clearly 
discussed. In particular, the red and blue curves in Fig.[T] 
are bounds on Tr from gravitino and axino production 
overlayed, assuming that they are co-dark matter candi- 
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lO""* 10-5 10-* 10-3 iQ-2 jo-1 lo" 10' 10^ 10^ 10-« 10-5 10-" 10-3 iQ-2 jo-1 lo" 10' 10^ 10' 

mi[GeV] [GeV] 



FIG. 2: Contours of fia/i^ = 0.11 for the KSVZ Q model, 
where cixinos couple minimally, for mj — 638 GeV. The 
{red, orange, yellow, green, blue} contours correspond to 
logio fa/GeV = {9, 10, 11, 12, 15}. 



dates and ma = TO3/2- This is a conservative assumption 
to make because as shown in Eq. (j?]), the axino acquires 
an irreducible contribution to its mass of order the grav- 
itino mass which cannot be eliminated without arbitrary 
fine-tuning. On the other hand, if for whatever reason the 
axino is heavier than the gravitino mass, then this only 
harshens axino bounds depicted in Fig.[Tl since the axino 
contour is translated to the left by a factor ma/m^/2. 

Fig. [T] illustrates a primary claim of this paper: the 
combined gravitino and axino problems entirely exclude 
the possibility of a high reheating temperature. The upper 
bound on Tr is about 3x 10^ GeV for fa = lO^^ GeV, and 
is even lower for lower fa- Thus, supersymmetry plus the 
axion solution to the strong CP problem are strongly at 
odds with theories of high Tr and high-scale leptogenesis, 
which typically require Tr > 10^ GeV. The bound on Tr 
is greatly relaxed for very large /a, as could occur with 
anthropic selection of a small axion misalignment angle. 



IV. GRAVITINO AND AXINO COSMOLOGY 

Finally, let us consider cosmology in theories with grav- 
itino LSP and axino NLSP, keeping careful track of the 
decays of axinos to gravitinos and cosmological limits on 
both axinos and saxions. The operator in Eq. Q gener- 
ates an axino-axion-gravitino coupling of the form 

Ga^&^d.a, (13) 

W3/2mpi 



FIG. 3: Contours of ilah'^ = 0.11 for the DFSZ [J] model, 
where axinos have additional couplings to squarks, for mg = 
638 GeV and = 500 GeV. The new decay q ^ a + q, 
results in a larger freeze-in region, where the contours are 
nearly vertical, as compared to Fig. (2] The contour colors 
correspond to the same values of fa as in Fig. [S] 

where G is the physical goldstino. Note that this operator 
respects the PQ symmetry because the axion is deriva- 
tively coupled. The above operator mediates the decay 
of the axino to the gravitino and axion with a lifetime 




This lifetime ranges over many orders of magnitude in 
the (m3/2, wa) plane, yielding a broad spectrum of cos- 
mological histories which we now explore. 



A. Three Scenarios for Dark Matter 

The cosmological history varies substantially in the 
(ma,TO3/2) plane, as depicted in Fig. [Hand Fig. [5] for 
fa = 10^ GeV and 10^^ GeV respectively The grey re- 
gion corresponds to ma < rn3/2, which is disfavored by 
the naturalness argument leading to Eq. ([3]), while the 
blue regions labeled by la, lb, II, and III are consistent 
with existing experiments. 

Region I of these figures, ra > 10^* sec, corresponds 
to a "Cosmologically-Stable Axino". In this regime, 
gravitinos and axinos are co-dark matter with mala + 
"^3/2^^3/2 = Pq/sq, where 

^ = 4 X 10-1" GeV (15) 
So 
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FIG. 4: Light blue regions in the (^3/2, n^a) plane have realis- 
tic dark matter cosmologies with a gravitino LSP and an axino 
NLSP for fa = 10^ GeV. Regions la and lb correspond dark 
matter comprised of both stable gravitinos and axinos, with 
gravitinos or axinos dominating the dark matter abundance, 
respectively. Region II is bounded by warm dark matter con- 
straints and limits on relativistic species. In Region III the 
axino mass is sufficiently large that the produced gravitinos 
are cold. The superpartner spectrum is taken to be identical 
to that of Fig. [T] 

is the observed dark matter relic abundance today. Here 
Ya and Y3/2 are the thermal axino and gravitino yields 
arising from scattering and decay processes. Provided 
that TOa,m3/2 > keV, both components are cold. Re- 
gion la corresponds to the dark matter energy density 
dominated by gravitinos, while Region lb corresponds to 
axino domination. 

Meanwhile, as the axino mass increases, the lifetime for 
the decay a ^ G + a decreases, so the axinos become cos- 
mologically unstable. Gravitinos from axino decay will 
then contribute to the total gravitino dark matter abun- 
dance, requiring 7713/2 (la + ^3/2) = Po/sq. 

If rria ^ ^13/2 the gravitinos produced in axino decay 
are highly relativistic, becoming non-relativistic only at 
a temperature 

Tnr^10-%Vx(^)'^'. (16) 

In Region III of Fig. |4] and Fig. [5l the axino mass 
nia ^ TeV, so that Tnr > keV and the gravitinos from 
axino decay are cold. There is little parameter space for 
this "Fast Decaying Axino" case, since axinos are are no- 
longer expected to be the NLSP if mg increases much 
above a TeV. 

For cosmologically unstable axinos, < TeV gives 
Tnr ^ keV so that the gravitinos from axino decay 



FIG. 5: Same as Fig. g] but with /„ = 10^^ GeV. Region 
II, corresponding to sub-dominant hot gravitino dark matter 
from axino decay, has increased. The superpartner spectrum 
is taken to be identical to that of Fig. [T] 

are warm and will over-deplete sub-Mpc structures un- 
less they constitute less than 3% of dark matter [s^]. 
The right-hand vertical contour that bounds Region II, 
the "Sub-Dominant Axino" case, corresponds to Ya < 
0.0313/2, placing an upper bound on 7713/2 which in- 
creases with fa as seen by comparing Fig. |1] and Fig. [5l 
The lower bound of Region II arises from the adverse 
effects of relativistic gravitinos arising from axino decay. 
As discussed in detail in [sOl, additional relativistic de- 
grees of freedom can adversely affect matter-radiation 
equality and are thus constrained by observations of the 
cosmic microwave background, galaxy clustering, and the 
Lyman-a forest. These bounds apply to axino lifetimes 
shorter than 10^"^ sec, and correspond roughly to an addi- 
tional effective number of neutrino species less than one. 
Numerically, this implies [soj 

™.n<3.4xl„-=G.Vx(-J^) {^).(U) 

where Ta is the temperature when the axino decays. Here 
Ta is simply computed by solving for the temperature 
at which H = l/ra, where ra is defined in Eq. (|14l) . 
Plugging in the value of Fa from scattering in Eq. (flT)) 
then yields a substantial limit on Tr given by 

For TTia ~ '713/2 the above bound on Tj^ is much more 
stringent than the one depicted in Fig. [T] from overclo- 
sure from gravitinos and axinos. However, by fixing Tr 
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to lie on the green contour of Fig. [1] such that gravitino 
dark matter from thermal production equals the mea- 
sured value today, then the above inequality for Tr can 
be reinterpreted as a lower bound on rria as a function 
of This constraint is what produces the slanted 

lower-right boundary of Region 11. 

In Fig. [B] contours of Q,^/2h^ = 0.11 have been plotted 
for several values of fa and a fixed superpartner spec- 
trum. These are contours of TO3/2 {Ya + ^5/2) = Po/so 
and therefore corresponds to gravitino dark matter in 
the "Sub-Dominant Axino" case. The contours become 
dashed when the values of 7/13/2 are excluded by the warm 
dark matter bound Ya < 0.031^3/2; this corresponds to 
the rightmost boundary of Region II in Fig. |4] and Fig. [5] 
Examining Fig. [S] it is clear that while the angled grav- 
itino scattering region is further excluded, the vertical 
gravitino freeze-in region is relatively unaffected. For 
fa = 10^^ GeV and 10^^ GeV the axino and gravitino 
yields in the bound Ya < 0.03y3/2 are both dominated 
by the scattering contributions. Thus the Tr depen- 
dence cancels and the resulting upper bound on 7713/2 
is dependent only on fa- As fa decreases the 7773/215 
contours move to left as seen in Fig. [21 Therefore the 
gravitino yield that is to be compared with the axino 
scattering yield becomes dominated by freeze-in, forcing 
Tr to smaller values. Since Yi"''^^ / Y^l'^'^^ cx ml^^TR/f^ 
the upper bound on relaxes whenever we cross over 
into the freeze-in dominated region. 

Fig. m and Fig. [5] depict cosmological restrictions on 
(™3/2i '^ai /a) which limit the range of collider signals 
resulting from the decay of the lightest observable-sector 
superpartner (LOSP) — that is, the lightest R-parity odd 
superpartner of a standard model particle. Assuming an 
axion model in which the LOSP couples directly to the 
axino (for instance a gluino LOSP in the KSVZ model) 
the ratio of LOSP decay rates to axinos and gravitinos is 
proportional to (7773/2 //a)^, then branching ratio to axi- 
nos (gravitinos) dominates at large (small) 7(73/2®. Hence 
in Region III at large 7713/2 the LOSP lifetime is fixed by 
fa and the dark matter is produced via axino production. 
At some intermediate values of 7773/2 the LOSP has siz- 
able branching ratios to both axinos and gravitinos, and 
if the axino mass is sufficiently heavy these modes can be 
distinguished by kinematics. Hence in parts of regions I 
and II it may be possible to measure (7773/2, rua, fa)- At a 
particular value of 7773/2, for instance around 7773/2 ~ 0.3 
MeV for squark masses of 500 GeV, the dark matter is 
produced by gravitino freeze- in. This value of 7/73/2 is 



Regardless of the choice of axion model any ultraviolet physics 
which generates the operator in Eq. JJjl will typically induce ax- 
ino/goldstino kinetic mixing via J t{A + A'')(X + X^), where 
e ~ fa/nipi. In components, eGa^^d^a^ is removed by shifting 
G — >■ G ea, which induces a coupling of the axino to the super- 
current of the form taj . Therefore, the LOSP typically decays to 
a with a branching fraction suppressed by a factor of relative 
to the branching fraction to G. 



10' - 
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FIG. 6: Gravitino dark matter in Region II, with contours 
of 7/73/2 (la -I- I3/2) = Po/sq. Here the {Red, Orange, Yellow, 
Green} contours correspond to logjQ /a/GeV = {9, 10, 11, 12}. 
The dotted portion of each contour shows values of 7/13/2 that 
are disallowed by hot dark matter, where hot gravitinos from 
axino decay constitute more than 3% of the total abundance. 
The transition from solid to dashed contour corresponds to 
crossing from Region II into the white region of Fig. U and 
Fig. [5] The superpartner spectrum is taken to be identical to 
that of Fig. [U 



sufficiently small that the LOSP decays dominantly to 
gravitinos and has a lifetime directly correlated with the 
freeze-in production mechanism 23]. For the charged 
slepton LOSP, decay signals to gravitinos and axinos have 
been studied in ^333 , examining the degree to which axino 
and gravitino modes can be distinguished. 



B. Cosmological Bounds from Saxions 

So far we have ignored the saxion component of the ax- 
ion supermultiplet. With communication between the su- 
persymmetry breaking and PQ sectors mediated at mass 
scale A, the saxion picks up a mass at the same order as 
the axino mass, mg ~ F/ K. This follows from very sim- 
ilar arguments to those of Section 2 — in particular, by 
inserting one power of Fi of Eq. ([9]) into the operator of 
Eq. ([7]), or more directly from the dimension 6 operator 
<i>j<i>.jA^A/A2 in the Kahler potential. In the effective 
theory beneath fa , Planck-scale dynamics induce the op- 
erator 

/ (f^e ^ '- ~ -7/73/2SS + h.c. + . . . ,(18) 

J 777p[ 2 
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so that the gravitino mass provides a lower bound for the 
saxion mass. 

The cosmological bounds imposed by saxions have 
been studied in detail in "s^. We find that under the 
assumption rus > 7713/2 these bounds on are typically 
less stringent than those coming from axino cosmology. 
In particular supersymmetric axion models must satisfy 
constraints on 

• Saxion overclosure, 

• Saxion decays infiuencing BBN, 

• Relativistic degrees of freedom. 

The saxion abundance receives a contribution from scat- 
tering in the early universe which is identical to the axino 
scattering yield yj^'^^.tt _ y^scatt^ g^g^^ from the 

supersymmetric interaction of Eq. ([l} . We assume a neg- 
ligible contribution to the saxion abundance from coher- 
ent oscillations. Unlike the axino and gravitino, saxion 
overclosure does not pose a threat because saxions decay 
rapidly; to photons through the electromagnetic analogue 
of Eq. (jlj and to axions through the kinetic term of the 
axion supermultiplet. 

Since the saxion will rapidly decay to photons one may 
worry that these photons might ruin the predictions of 
BBN unless the saxion is made heavy enough. However, 
this is not the case, since decays to photons are sup- 
pressed by a loop factor and so will always be subdomi- 
nant to decays to axions 

r ^ aa) _ 64712 ^^^^ 
r (s — > 77) ' 

where a is the fine structure constant. While the branch- 
ing fraction of the saxion to axions can in principle be 
small, this requires a delicate and unnatural cancellation 
in the underlying theory. 

Finally, let us consider the issue of relativistic degrees 
of freedom. In complete analogy with the axino, the sax- 
ion will decay and produce relativistic energy in axions 
which is subject to the bound in Eq. (|17l) taken from [30] 
except with maYa and Ta replaced with mgYs and T,. 
Plugging in for Yg, the saxion yield from thermal scat- 
tering (which is essentially equal to the axino yield from 
thermal scattering) , and plugging in for Tg , the tempera- 
ture at the time of saxion decay where ~ m^/6ATrf^, 
one acquires a bound on Tji given by 

7-.<^xl»'GoVx(^)(^)'^ 

which is generally weaker than the Tr bound from over- 
closure depicted in Fig. [1] 



LSP and axino NLSP. In this case the axino plays an im- 
portant role in determining viable cosmological histories. 
In particular, the overclosure constraint alone provides a 
very powerful limit on the reheat temperature, as shown 
in Figure 6, that can be approximated by 

^-<^^io'^^^(to4;v)- 

We identify three very different phases of combined grav- 
itino and axino cosmology which we label according to 
the nature of the axino: 

(I) "Cosmologically-Stable Axino" is the case 
where axinos and gravitinos are co-dark matter. 

(II) "Sub-Dominant Axino" production leads to 
gravitino dark matter from decays, with an upper 
bound on 7713/2 that depends on fa- 

(III) "Fast-Decaying Axino" has a TeV scale mass 
and thus decays sufficiently early that gravitino 
dark matter is cold. 

The regions of gravitino and axino masses that yield 
these three cosmologies are shown in Fig.|3]and Fig.[5]for 
fa — 10^ GeV and 10^^ GeV. In all three cosmologies, 
dark matter may be dominantly produced by gravitino 
freeze-in. 

The required ranges of Tr for cases (II) and (III) are 
shown in Fig. [5] for several values of fa- They give an 
upper bound on Tr that is much lower than for the 
case of gravitino dark matter without axinos, and is also 
lower than the bound of Eq. ([^0]) . Thus the presence of 
the axinos increases the likelihood that the production 
mechanism for gravitino dark matter is freeze- in from de- 
cays, strengthening the possibility that LHC will provide 
strong evidence for gravitino dark matter by measuring 
the LOSP lifetime [2^. Finally, let us note that our con- 
clusions are unaltered by cosmological considerations of 
the saxion for any value of its mass greater than 7773/2- 
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Appendix A: A Simple Theory 



CONCLUSIONS 



Consider a minimal theory of PQ and supersymmetry 
breaking defined by 



If supersymmetry and the QCD axion coexist, then a £ f d^O [ d^O h (A 

large range of parameter space will include a gravitino ~ / ^ ^ J ^+ -C-j ( 1) 
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where the Kahler potential and and superpotential are 
defined by 



K ^ G{^^<i>) + H{X\X) + 



A 



W = FX. 



h.c. (A2) 
(A3) 



Here $ and X are the PQ and supersymmetry breaking 
fields, respectively, and G and H are as of yet unspec- 
ified real functions. Since $ is PQ charged, it cannot 
be present in the superpotential. Also, note that G is 
purely a function of so that the PQ symmetry is 
left unbroken. We interpret the 1 /A suppressed operator 
coupling $ and X as something generated by unspecified 
high-scale dynamics. 

Demanding that G and H have a form such that (<i>) = 
fa and {X) = at the minimum of the scalar potential 
implies that 



(G(i)) = fAG'-'^) + f^{G<''^) (A4) 



J a 



A3((G(i)) + /2(G(2)))2 



where the superscripts denote derivatives with respect to 
the function arguments. Plugging these expressions into 
the action, canonically normalizing the fermion kinetic 
terms, and then diagonalizing the fermion mass matrix, 
one discovers that one linear combination of fermions 
from $ and X is massless, as is expected of the goldstino. 
Meanwhile, the orthogonal axino component acquires a 
mass 



F 



m~a 



2A (G(i)) -^/2(G(2)) (iJ(i^i) 



(A5) 



where the ellipses denote terms higher order in A. Hence, 
we confirm the result of the simple operator argument in 
Eq. ©. 
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